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By Haecheon Choi and Roger Temam 1 


In this article, we make a few remarks on the role that attractors and inertial 
manifolds play in fluid mechanics problems. We then describe the role of incremental 
unknowns for approximating attractors and inertial manifolds when finite difference 
multigrid discretizations are used. The relation with direct numerical simulation 
and large eddy simulation is also mentioned. 

1. Motivations and objectives 

At this time, there are two different methods for mathematically describing a 
turbulent flow when the permanent regime is established. The first one, more tra- 
ditional, is related to the concept of ensemble averages and the idea that turbulent 
flows are statistically well defined and reproducible. The flow is described by a 
measure in the infinite dimensional function space, the statistics being, of course, 
that associated to the measure. Probabilities and statistical tools are essential in 
this approach. 

Another more recent approach to turbulence stems from dynamical system the- 
ory in which the flow is described by an attractor which can be a complicated 
(fractal) set (see e.g. Constantin, Foias h Temam 1985). Furthermore, this at- 
tractor is expected to have in general a large dimension although some interesting 
turbulent flows with a low dimensional attractor have been produced in special 
cases. However, this low dimensional behavior is not likely to appear in general for 
flows of industrial interest (see Keefe, Moin & Kim 1992). Before making further 
remarks on the attractor point of view, let us observe that the statistical and at- 
tractor approaches are not inconsistent. A measure can be defined on the attractor. 
Furthermore, it was one of the main results of Constantin, Foias & Temam (1985, 
1988) to show that certain segments of the Kolmogorov and Kraichnan theories of 
turbulence can be rigorously derived from the Navier-Stokes equations using the 
attractor point of view. 

It may not seem practically useful to state that a turbulent flow is represented 
by a fractal attractor of large dimension. Fortunately, some useful information 
can be derived from the study of the attractor which sheds some new light on the 
numerical simulation of turbulence. The first information is that the attractor can 
be embedded approximately, or perhaps exactly, in smooth manifolds called inertial 
manifolds (Foias, Sell & Temam 1985, 1988; Foias, Manley & Temam 1987; see also 
Temam 1991a). These manifolds yield a slaving of the high frequency component of 
the flow by its low frequency component and reduce the number of modes needed 
to monitor a flow. 
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Incremental unknowns have been introduced in Temam (1990) as a means to 
approximate attractors and inertial manifolds in the context of finite difference 
multigrid approximations. Indeed, when spectral discretization is used, the separa- 
tion between the small and large scale components of the flow occurs naturally in 
the spectral space. When finite differences are used, all mesh points play the same 
role however large is the number of mesh points, and there is no obvious way to 
distinguish between the small scales which carry little energy and the large scales 
which carry most of the energy. Incremental unknowns of various types have been 
introduced which produce the separation of scale which is needed (Chen & Temam 
1993). The utilization of incremental unknowns for large scale computations of tur- 
bulent flows, for direct numerical simulation or large eddy simulation, remains to be 
done. As a first step towards this important task, we describe hereafter in section 
2 the utilization of incremental unknowns for the solution of the Burgers equation 
with stochastic forces. Indeed, these equations are known to be a good model for 
turbulent flows (see e.g. Chambers ct ah 1988; Choi et ah 1992). 

These computations support two important observations: 

The small scale component of the flow as defined in the context of incremental 

unknowns (IU) is indeed small. Its instantaneous variation is very fast while its 
averaged variation is very slow. 

The IU method is numerically efficient. It produces an improvement of the 

CFL stability condition which has not been analytically fully explained, although 
there are already some partial theoretical justifications (Temam 1990, 1991b). Con- 
sequently, the same accuracy can be recovered with less calculations, the gain in 
computing time in the cases that we consider being a factor of 7. 


2. Accomplishments 

In this section, we describe the procedure for applying the IU method to the 
stochastic Burgers equation. This equation contains nonlinear convection and dif- 
fusion terms, and its solution exhibits a chaotic nature; these qualities make it a 
natural model for the more complicated Navier-Stokes equations. Section 2.1 de- 
scribes the governing equation and boundary conditions. The mathematical and 
numerical procedures of applying the IU method to the stochastic Burgers equation 
are described in section 2.2. Numerical results are presented in section 2.3. Ap- 
plication of the IU method to the Navier-Stokes equations is considered in section 
2.4. 

2.1 The Burgers equation with random forcing 
Consider the randomly forced non-dimensionalized Burgers equation with no-slip 
boundary conditions 


du d u 2 
~dt + foT 


1 d 2 u 
Re dx 2 


+ x(M), 


0 < x < 1, 


u(x = 0 ,t) = u(x = 1,<) = 0, 
u(x, t = 0) = u 0 , 


( 2 . 1 ) 
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where u is the velocity, \ is the random forcing, and Re is the Reynolds number. 
The initial data u Q is an instantaneous solution of the Burgers equation with random 
forcing \ rather than an arbitrary function. The forcing function \ is a white noise 
random process in x with zero mean (see Chambers et al. 1988; Choi et al. 1992): 

< X >1= 0* < X 2 >*= 1- 

Here < • > x denotes the average value over space. In the absence of forcing (x “ 0), 
the solutions of equation (2.1) decay to zero from any bounded initial data. 

2.2 Description of the IU method for the Burgers solution 

The discretized equations for the stochastic Burgers equation with no-slip bound- 
ary conditions are obtained using a semi-implicit method in time (Adams-Bashforth 
for the convection term and Crank- Nicholson for the diffusion term) and a second- 
order centered difference method in space: 

au" - /?u " +1 - = < 7 ,, * = 1 , • • *, 2 / — 1 , ( 2 . 2 ) 


«o = u 2 i = 0, 


where 


a = 1 + 
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2ReAx 2 
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o/ n — l 2 — 1 2 \ / 71 — 22 71— 2 2 \ , A J 

3( u t+l ~ u i- 1 ) ( u i+l u i-l i 


(2.3) 


where At is the computational time step and Ax = 1/(2 1) is the grid spacing. 

The incremental unknowns for the present problem consist of the numbers <f> 2 i 
and V>2*+i (see figure 1): 


(f>2i — ^2i? i — 9, * * ■, JT, 

02i+l = ^2i + l “ 2^ 2 * * = 0, •••,/- 1, 


(2.4) 


where <f)$ and (f> 2 i correspond to the velocities at the boundary points. Thus ip 2 %+\ 
is the increment of u to the average of the values at the neighboring points, 2 i and 
2 i + 2; hence by Taylor’s formula, V^i+i is small, of order Ax 2 . 

At points 2i + 1, equation (2.2) becomes 


au 2t+l — /^ w 2t+2 0 u 2i ~ 02i+l- 


132 


H. Choi & R. Temam 



FIGURE 1. Incremental unknowns. 


Using equation (2.4) gives 

«02i+l + (f " (^" I + ^" 1 + 2 ) = (2-5) 

where i = 0, •••,/ — 1. Similarly, at points 2i, equation (2.2) becomes 


aw 2i — ^ W 2i+1 /^ w 2i-l — ^2i* 

Substituting equations (2.4) and (2.5) into the above equation gives 

^ (« 2 - 2/3 2 ) <^2i - ^+",+2 “ ~^2i-2 = 92i + ^ (fl^2t+i + <721-1 ) , (2.6) 

where i = 1,. • •,/ — 1. The system consisting of equation (2.6) is similar to the 
system consisting of equation (2.2), but involves half as many unknowns. 

We can, of course, repeat the procedure. If we start with Ax = 1/(2 '/), then 
after l steps, we reduce the initial system involving 2 l I unknowns to a similar one 
involving I unknowns. Following the same procedure described above, one can get 
the discretized equations for multi-steps 1: 


m jn,m nm ifi.m /U,m 771 - 1 t <v~ 7 n r __ 1 

& Y2i 0*2 1 + 2 P 02i— 2 — #2i 5 Z A A 7 


(2.7) 


; n . m 

^2i+l = 


a' 


^ + ,-(%--r)(0 2 r+« 2 ) 


, i = 0, • • — 1, (2.8) 


a " = ( am |2 — 2/? m ~ l2 ) , 


(i m = 


F 


i-l2 


CV 


m — 1 


(2.9) 


J17I — 1 
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where m = 1, ■ • / denotes the level of grids and m = 1 and l correspond to the 

finest and coarsest grids, respectively; 

a 0 = a, /3° = /? and g? = i = 1, • * *, 2 l I - 1 (equation (2.3)). 

Numerical procedure 

The numerical procedure of applying the IU method to the stochastic Burgers 
equation can be written as follows: 

Step 1: Start from an initial velocity field u 0 , u”"* 1 = 

Step 2: Get a m , /3 m , and g from equation (2.9) for m = 1, ■ ■ L 
Step 3: Solve equation (2.7) for the coarsest grid (m = /) to obtain tp”?- 
Step 4: Obtain V^i+i f rom equation (2.8) with 
Step 5: Obtain u n,m ; u™'™ = <p 
Step 6: Obtain 

Step 7: Repeat Steps 4 and 6 until u n>l is obtained. 

2.3 Numerical results using the IU method 
In this section, we apply the numerical procedure described in the previous section 
to the Burgers equation. A uniform computational mesh of 2049 points is used in 
x (Ax = 1/2048) and l — 4 (number of the grid levels); therefore, the coarsest grid 
has 257 points. 

2.3.1 Properties of the incremental unknowns 

An initial velocity field u 0 is obtained for Re = 1500 and A t r = 0.01, where A t r is 
the time scale of the random forcing (for more details see Choi et al. 1992). Figure 
2 (a) shows the root-mean-square values of <f> and tp m as a function of where 
m = 1,2,3; m = 1,2, and 3 correspond to the grid points of 2049, 1025, and 513, 
respectively. It can be seen that the 0 m, s are several orders of magnitude smaller 
than <f>. Also, note that the incremental unknowns tp at the coarser grid level have 
larger magnitude as compared to those at the denser grid level. 

The root-mean-square values of d<p/dt and dtp™ /dt as a function of t are shown 
in figure 2 ( b ). Contrary to the results of <j> and tp™ (figure 2 (a)), the magnitude of 
dtp™ /dt is comparable to that of d<p/dt, indicating that dtp™ fdt cannot be neglected 
pointwise as compared to d<f>jdt when there is a stochastic motion in flow. In 
case there is no small-scale motion in flow, however, the magnitude of dip™ /dt is 
several orders of magnitude smaller than that of d<f>/dt : we have tested an initial 
velocity field, sin ( 27 rx), without random forcing, and the result showed this is indeed 
the case. When a small-scale motion exists in the flow, we expect the derivatives 
dtp™ /dt to be small in average as evidenced in figure 3. This point will be discussed 
elsewhere. 

Figure 4 shows dtp 1 /dt at x = 0.5 as a function of t and dtp 1 fdt at t = 2 as a 
function of x; dtp* /dt(t) shows an intermittent behavior while a chaotic behavior is 
shown in dtp* /dt(x). The power spectra of these dtp* fdt's show that the incremental 
unknowns have nearly same power at all wavenumbers (a white noise) while they 
have higher power at the highest frequency. 
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(b) 


X 

FIGURE 4. Instantaneous as a function of t or x. 


2.3.2 Efficiency of the IU method 

In this section, we investigate the effect of the IU method on the numerical 
stability. A classical semi-implicit method has a limit of the computational time step 
due to the explicit treatment of the nonlinear term. The maximum C FL (Courant- 
Friedrichs-Lewy) number is restricted to one with the present semi-implicit method 
(equation (2.2)). Hence, the maximum time step without the IU method is restricted 
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Case 

l = 1 

1 = 2 

1 = 3 

1 = 4 

At ma x 

0.0051 

0.017 

0.054 

0.054 

Total CPU (sec) 

0.232 

0.088 

0.032 

0.034 


Table 1. Maximum computational time step and total CPU time using the IU 
method. 


( 2 . 10 ) 

Maximum computational time steps At max using the IU method have been ob- 
tained for the initial value problem (2.1): u Q — sin ( 2nx ), Re = 100, and \ = 0 ( no 
random forcing). Table 1 shows the maximum computational time step and total 
CPU time using the IU method, where l denotes the level of the coarsest grid. Note 
that the IU method is not used when / = 1. A maximum computational time step 
is clearly increased when / > 1, indicating that the numerical stability is enhanced 
using the IU method. However, there seems to exist a limit of level l of increasing 
At ma x\ for the present problem, l = 3 (see table 1). 

The CPU time to advance one time step with the IU method is larger than 
that without the IU method due to the calculation procedure of and g m , 

although the matrix size to be inverted is reduced from 2 l I to I. Total CPU times 
to reach t = 1 are shown in table 1. Computational cost is significantly reduced 
by a factor of 7 using the IU method. The mathematical analysis of the stability 
enhancement using the IU method for the Burgers equation is very complex and will 
be investigated in the future. Theoretical indications that the IU method produces 
stability enhancement appear in Temam (1990, 1991). 

2.4 Application of the IU method to the unsteady Navier-Stokes equations 

The IU method can be applied to the unsteady three-dimensional Navier-Stokes 
equations as follows: Firstly, the first step of the fractional-step method (Kim &; 
Moin 1985) provides a second-order-accurate approximation of the three-dimensional 
Navier-Stokes equations with pressure term excluded. Secondly, an approximate 
factorization technique (Beam & Warming 1978) splits three-dimensional equations 
to three one-dimensional equations. Finally, the IU method as presented above is 
applied to each one-dimensional equation. This case will be addressed in the future. 
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